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Abstract 

A systematic study of the factorization of long-range azimuthal two-particle corre¬ 
lations into a product of single-particle anisotropies is presented as a function of pj 
and ?/ of both particles, and as a function of the particle multiplicity in PbPb and 
pPb collisions. The data were taken with the CMS detector for PbPb collisions at 
^/s NN = 2.76 TeV and pPb collisions at = 5.02 TeV, covering a very wide range 
of multiplicity. Factorization is observed to be broken as a function of both particle pj 
and //. When measured with particles of different pj, the magnitude of the factoriza¬ 
tion breakdown for the second Fourier harmonic reaches 20% for very central PbPb 
collisions but decreases rapidly as the multiplicity decreases. The data are consistent 
with viscous hydrodynamic predictions, which suggest that the effect of factorization 
breaking is mainly sensitive to the initial-state conditions rather than to the trans¬ 
port properties (e.g., shear viscosity) of the medium. The factorization breakdown is 
also computed with particles of different;/. The effect is found to be weakest for mid¬ 
central PbPb events but becomes larger for more central or peripheral PbPb collisions, 
and also for very high-multiplicity pPb collisions. The //-dependent factorization data 
provide new insights to the longitudinal evolution of the medium formed in heavy 
ion collisions. 
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1 Introduction 


The goal of experiments with heavy ion collisions at ultra-relativistic energies is to study nu¬ 
clear matter under extreme conditions. By studying the azimuthal anisotropy of emitted par¬ 
ticles in such collisions, experiments at the Relativistic Heavy Ion Collider at BNL (RHIC) in¬ 
dicated that a strongly-coupled hot and dense medium is created, which exhibits a strong col¬ 
lective flow behavior IM. At the significantly higher collision energies achieved at the Large 
Hadron Collider (LHC), the collective phenomena of this quark gluon plasma have also been 
studied in great detail 1)511131. 

The collective expansion of the hot medium in heavy ion collisions can be described by hy¬ 
drodynamic flow models. Motivated by such models, the azimuthal distribution of emitted 
particles can be characterized by the Fourier components of the hadron yield distribution in 
azimuthal angle (<p) ffl4hT6ll , 

^ oc 1 +2^p„cos[n((/> - Y„)]. (1) 


Here, the Fourier coefficients, v n , characterize the strength of the anisotropic flow, while the 
azimuthal flow orientation is represented by the corresponding "event plane" angle, Y„, the 
direction of maximum final-state particle density. The event plane angles are related to the 
event-by-event spacial distribution of the participating nucleons in the initial overlap region. 
The most widely studied and typically also strongest form of anisotropic flow is the second 
Fourier component, Vi, called "elliptic flow". The elliptic flow event plane, Y 2 , is correlated 
with the "participant plane" given by the beam direction and the shorter axis of the approxi¬ 
mately elliptical nucleon overlap region. Because of event-by-event fluctuations, higher-order 
deformations or eccentricities of the initial geometry can also be induced, which lead to higher- 
order Fourier harmonics ( v n , n > 3) in the final state with respect to their corresponding event 
plane angles, Y„ H7l . Studies of azimuthal anisotropy harmonics provide important informa¬ 
tion on the fundamental transport properties of the medium, e.g., the ratio of shear viscosity to 
entropy density, ;//s Iil8] - t20l . 

A commonly used experimental method to determine the single-particle azimuthal anisotropy 
harmonics, v n , is the measurement of two-particle azimuthal correlations lil4tiT6l : 21 |. The az¬ 
imuthal distribution of particle pairs as a function of their relative azimuthal angle Acp can also 
be characterized by its Fourier components. 


dN pair 
dA (p 


oc1 + 2£Ka cos(hA^). 

n 


( 2 ) 


If the dominant source of final-state particle correlations is collective flow, the two-particle 
Fourier coefficients, V n ,\, are commonly expected to follow the factorization relation: 

VnA = K v b nr (3) 

where v" and v b n represent the single-particle anisotropy harmonics for a pair of particles (a 
and b) in the event. The particle pairs can be selected from the same or different transverse 
momentum (pj) and pseudorapidity (;/) ranges. Here, a key assumption is that the event plane 
angle Y„ in Eq. ([T]) is a global phase angle for all particles of the entire event, which is canceled 
when taking the azimuthal angle difference between two particles. As a result, the flow-driven 
A(p distribution in Eq. (|2j has no dependence on Y„. The most common approach to obtain the 
single-particle v n in the two-particle method is to fix one particle in a wide pj (rj) region and 
measure V n & by only varying pj (tj) of the other particle to determine v„ as a function of pj {rj). 
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1 Introduction 


However, a significant breakdown of the factorization assumption, up to about 20%, was re¬ 
cently observed for pairs of particles, separated by more than 2 units in //, from different pj 
ranges in ultra-central (0-0.2% centrality) PbPb collisions Itl3l . The centrality in heavy ion 
collisions is defined as a fraction of the total inelastic PbPb cross section, with 0% denoting 
the most central collisions. While nonflow correlations (such as back-to-back jets) have been 
speculated to possibly account for this effect, contributions of those short-range correlations to 
the collective anisotropy are less dominant in high-multiplicity events as the total number of 
particles increases ||22| . It was then realized that in hydrodynamic models the assumption of 
factorization does not hold in general because of fluctuations in the initial overlap region of 
two nuclei |23:.:24J- In each event, due to local perturbations in the energy density distribution 
generating a pressure gradient that drives particles in random directions with differing boosts, 
the resulting event plane angles found with final-state particles from different pj ranges may 
fluctuate with respect to each other (although still correlated with the initial participant plane). 
This effect of initial-state fluctuations thus breaks the factorization relation of Eq. (|3jl, which as¬ 
sumes a unique event plane angle for all particles in an event. As a result, the precise meaning 
of previous single-particle v„ results should be reinterpreted as being with respect to the event 
plane determined with particles over a specific, usually wide, pj range. Quantitative studies 
of the factorization breakdown effect as a function of pj could place stringent constraints on 
the spatial scale (or granularity) of the fluctuations in the initial state of heavy ion collisions, 
especially along the radial direction Il25r(27l . 

The recent observation of long-range near-side (A (p ~ 0) two-particle correlations in pp |28| 
and pPb ll29l - f3H collisions raised the question of whether hydrodynamic flow is developed 
also in these small collision systems. The extracted v n harmonics in pPb collisions have been 
studied in detail as a function of pj and event multiplicity [22, j32]. The initial-state geometry 
of a pPb collision is expected to be entirely driven by fluctuations. If the observed long-range 
correlations in such collisions indeed originate from hydrodynamic flow, the effect of factoriza¬ 
tion breakdown should also be observed in the data and described by hydrodynamic models. 
Since the initial-state geometries of both high-multiplicity pPb and ultra-central PbPb colli¬ 
sions are dominated by fluctuations, it is of great interest to investigate whether the magnitude 
of factorization breakdown is similar in these two systems. 

Furthermore, the factorization breakdown in // is sensitive to event plane fluctuations at dif¬ 
ferent // |23| . This phenomenon has been investigated in hydrodynamic and parton transport 
models Il33if36ll . The observation and study of this effect will provide new insights into the dy¬ 
namics of longitudinal expansion of the hot quark and gluon medium, and serves as an ideal 
test ground of three-dimensional hydrodynamic models. 

This paper presents a comprehensive investigation of the factorization breakdown effect in 
two-particle azimuthal Fourier harmonics in PbPb (pPb) collisions at = 2.76 (5.02) TeV, 
to search for evidence of pj- and //-dependent event plane fluctuations. The Fourier harmonics 
of two-particle azimuthal correlations are extracted for pairs with |A//| > 2 as a function of pj 
and // of both particles in a pair. The results are presented over a wide range of centrality or 
event multiplicity classes, and are compared with hydrodynamic models in PbPb and pPb col¬ 
lisions. As the pj- and //-dependent aspects of factorization breakdown probe system dynamics 
in the transverse and longitudinal directions, respectively, an assumption is made that the de¬ 
pendence on each variable can be studied independently by averaging over the other, and two 
different analysis techniques are applied. These two aspects of the analysis are described in 
Sections [4] and [5] separately, including the analysis procedures and results. 
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2 Experimental setup and data sample 

A comprehensive description of the Compact Muon Solenoid (CMS) detector at the CERN 
LHC, together with a definition of the coordinate system used and the relevant kinematic vari¬ 
ables, can be found in Ref. 1371 . The main detector sub-system used in this paper is the tracker, 
located in a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 
3.8 T. The tracker consists of 1440 silicon pixel and 15148 silicon strip detector modules, cover¬ 
ing the pseudorapidity range \ rj\ < 2.5. For hadrons with pj ~ 1 GeV/c and \rj\ ~ 0, the impact 
parameter resolution is approximately 100 pm and the pj resolution is 0.8%. 

The electromagnetic calorimeter (ECAL) and the hadron calorimeter (HCAL) are also located 
inside the solenoid. The ECAL consists of 75 848 lead tungstate crystals, arranged in a quasi- 
projective geometry and distributed in a barrel region (\tj\ < 1.48) and two endcaps that extend 
to |;/| = 3.0. The HCAL barrel and endcaps are sampling calorimeters composed of brass and 
scintillator plates, covering \rj\ < 3.0. In addition, CMS has an extensive forward calorimetry, in 
particular two steel/quartz-fiber Cherenkov hadronic forward (HF) calorimeters, which cover 
the pseudorapidity range 2.9 < \rj\ < 5.2. The HF calorimeters are segmented into towers, each 
of which is a two-dimensional cell with a granularity of 0.5 in ?/ and 0.349 radians in cp. A set 
of scintillator tiles, the beam scintillator counters (BSC), are mounted on the inner side of the 
HF calorimeters and are used for triggering and beam-halo rejection. The BSCs cover the range 
3.23 < |j/ 1 < 4.65. The detailed Monte Carlo (MC) simulation of the CMS detector response is 
based on Geant 4 Il38l . 

The data sample used in this analysis was collected with the CMS detector during the LHC 
PbPb run in 2011 and pPb run in 2013. The total integrated luminosity of the data sets is about 
159 /fb 1 for PbPb, and 35 nb 1 for pPb. During the pPb run, the beam energies were 4TeV 
for protons and 1.58 TeV per nucleon for lead nuclei, resulting in a center-of-mass energy per 
nucleon pair of 5.02 TeV. As a result of the energy difference between the colliding beams, 
the nucleon-nucleon center-of-mass in the pPb collisions is not at rest in the laboratory frame. 
Massless particles emitted at r/ an = 0 in the nucleon-nucleon center-of-mass frame will be 
detected at;/ = —0.465 or 0.465 (clockwise or counterclockwise proton beam) in the laboratory 
frame. 

3 Selection of events and tracks 

Online triggers, offline event selections, track reconstruction and selections are identical to 
those used in previous analyses of PbPb and pPb data llT3l |22| and are briefly outlined in the 
following sections. 

3.1 PbPb data 

Minimum bias PbPb events were selected using coincident trigger signals from both ends of 
the detector in either BSCs or the HF calorimeters. Events due to detector noise, cosmic rays, 
out-of-time triggers, and beam backgrounds were suppressed by requiring a coincidence of the 
minimum bias trigger with bunches colliding in the interaction region. The trigger has an ef¬ 
ficiency of (97 ± 3)% for hadronic inelastic PbPb collisions. Because of hardware limits on the 
data acquisition rate, only a small fraction (2%) of all minimum bias events were recorded (i.e., 
the trigger is "prescaled"). To enhance the event sample for very central PbPb collisions, a ded¬ 
icated online trigger was implemented by simultaneously requiring the HF transverse energy 
(Ej) sum to be greater than 3260 GeV and the pixel cluster multiplicity to be greater than 51400 
(which approximately corresponds to 9500 charged particles over 5 units of pseudorapidity). 
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3 Selection of events and tracks 


The selected events correspond to the 0.2% most central PbPb collisions. Other standard PbPb 
centrality classes presented in this paper are determined based on the total energy deposited 
in the HF calorimeters IfTTH . The inefficiencies of the minimum bias trigger and event selection 
for very peripheral events are properly taken into account. 

To further reduce the background from single-beam interactions (e.g., beam-gas and beam- 
halo), cosmic muons, and ultra peripheral collisions that lead to the electromagnetic breakup 
of one or both Pb nuclei Il39l , offline PbPb event selection criteria lfTT| are applied by requiring 
energy deposits in at least three towers in each of the HF calorimeters, with at least 3 GeV of 
energy in each tower, and the presence of a reconstructed primary vertex containing at least 
two tracks. The reconstructed primary vertex is required to be located within ±15 cm of the 
average interaction region along the beam axis and within a radius of 0.02 cm in the transverse 
plane. Following the procedure developed in Ref. Itl3l . events with large signals in both Zero 
Degree Calorimeter (ZDC) and HF are identified as having at least one additional interaction, 
or pileup events, and thus rejected (about 0.1% of all events). 

The reconstruction of the primary event vertex and of the trajectories of charged particles in 
PbPb collisions are based on signals in the silicon pixel and strip detectors and described in 
detail in Ref. Bill . From studies based on PbPb events simulated using hydjet vl.8 EOl . the 
combined geometrical acceptance and reconstruction efficiency of the primary tracks is about 
70% at pj ~ 1 GeV/c and \rj\ < 1.0 for the most central 0-5% PbPb events, but drops to about 
50% for pj ~ 0.3 GeV/c. The fraction of misidentified tracks is kept at the level of <5% over 
most of the pj (>0.5 GeV/c) and \rj\ (<1.6) ranges. It increases to about 20% for very low pj 
(<0.5GeV/c) particles in the forward (\t]\ > 2) region. 

3.2 pPb data 

Minimum bias pPb events were selected by requiring that at least one track with pj > 0.4 GeV/c 
is found in the pixel tracker in coincidence with a pPb bunch crossing. About 0.1% of all min¬ 
imum bias pPb events were recorded. In order to select high-multiplicity pPb collisions, a 
dedicated high-multiplicity trigger was implemented using the CMS level-1 (LI) and high- 
level trigger (HLT) systems. At LI, the total transverse energy measured using both ECAL and 
HCAL is required to be greater than a given threshold (20 or 40 GeV). Online track reconstruc¬ 
tion for the HLT was based on the three layers of pixel detectors, and required a track origin 
within a cylindrical region, centered at the average interaction point of two beams, of length 
30 cm along the beam and radius 0.2 cm perpendicular to the beam. For each event, the vertex 
reconstructed with the highest number of pixel tracks was selected. The number of pixel tracks 

(jyonUne) with |^| 

< 2.4, pj > 0.4 GeV/c, and a distance of closest approach of 0.4 cm or less to 
this vertex, was determined for each event. 

Offline selections similar to those used for the PbPb data sample are applied to reject non- 
hadronic pPb interactions. A coincidence of at least one HF calorimeter tower with more than 
3 GeV of total energy in each of the HF detectors is required. Events are also required to contain 
at least one reconstructed primary vertex within 15 cm of the nominal interaction point along 
the beam axis and within 0.15 cm transverse to the beam trajectory. At least two reconstructed 
tracks are required to be associated with the primary vertex. Beam-related background is sup¬ 
pressed by rejecting events for which less than 25% of all reconstructed tracks are of sufficiently 
good quality to be tracks selected for physics analysis, as will be discussed later in this section. 
Among those pPb interactions simulated with the EPOS Pill and HIJING B42l event genera¬ 
tors that have at least one primary particle with total energy E > 3 GeV in both // ranges of 
—5 < rj < —3 and 3 < // < 5, the above criteria are found to select 97-98% of the events. Pileup 
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events are removed based on the number of tracks associated with each vertex in a bunch cross¬ 
ing and the distance between different vertices Il22i . A purity of 99.8% for single pPb collision 
events is achieved for the highest multiplicity pPb interactions studied in this paper. 

For the pPb analysis, the standard track reconstruction as in pp collisions is applied. The CMS 
high-purity tracks (as defined in Ref. Il43l ) are used. Additionally, a reconstructed track is only 
considered as a primary-track candidate if the significance of the separation along the beam 
axis (z) between the track and primary vertex, d z /cr(d z ), and the significance of the impact 
parameter relative to the primary vertex transverse to the beam, dj/cr(dj), are each less than 
3. The relative uncertainty in the transverse momentum measurement, cr(pj) /pj, is required 
to be less than 10%. To ensure high tracking efficiency and to reduce the rate of misidentified 
tracks, only tracks within | t/ \ < 2.4 and with pj > 0.3 GeV/c are used in the analysis. 

The entire pPb data set is divided into classes of reconstructed track multiplicity, N°^ me , where 
primary tracks with \rj\ < 2.4 and pj > 0.4GeV/c are counted. The multiplicity classification 
in this analysis is identical to that used in Ref. ||22j, where more details are provided. The 
more central (0-50%) PbPb data, including ultra-central triggered events, are analyzed with 
a standard reconstruction algorithm used in heavy ion collisions, as described in Section |3.1| 
In order to compare the pPb and PbPb systems at the same collision multiplicity, peripheral 
PbPb events for 50-100% centrality are reprocessed using the same event selections and track 
reconstruction as for the pPb analysis. 


4 Transverse momentum dependence of factorization breakdown 

4.1 Analysis technique 

The pj-dependent factorization breaking effect is investigated using the same analysis tech¬ 
nique of two-particle azimuthal correlations as that applied in Ref. Ifl3l . For simplicity, a pair 
of two charged tracks are labeled as particle a and b (equivalent to the trigger and associated 
particles used in previous publications). They are selected from the same or different p\ and 
Pj ranges within \tj a,b \ < 2.4. The two-particle Fourier coefficients, V n/ \, are calculated as the 
average value of cos(nA^)) over all particle pairs, which fulfill the requirement of |A//j >2 (to 
avoid the short-range correlations from jets and resonance decays): 

V„A = ((cos(nA0)))s - ((cos(nA0))) B , (4) 

in given ranges of pj and p\. Here, (()) denotes averaging over all particle pairs in each event 
and over all the events. The subscript S corresponds to the average over pairs taken from the 
same event, while B represents the mixing of particles from two randomly-selected events in 
the same 2 cm wide range of the primary vertex position in the z direction and from the same 
centrality (track multiplicity) class. The ((cos(nA(/)))) g term, which is typically two orders of 
magnitude smaller than the corresponding S term, is subtracted to account for the effects of 
detector non-uniformity. This analysis is equivalent to those in Refs. [10, 22} |44j @5], where 
the two-particle azimuthal correlation function is first constructed and then fit with a Fourier 
series. The advantage of the present approach is that the extracted Fourier harmonics will not 
be affected by the finite bin widths of the histogram in At/ and A <p of the two-particle correlation 
function, which is relevant for higher-order Fourier harmonics. 
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4 Transverse momentum dependence of factorization breakdown 


With the V n &(pj, p\) values as a function of Pj and p\, the factorization ratio. 


r n{Pj, Pt) 


V»a(Pt/ Pt) 

\JVnA(Pj, P^VnAiPj, v\) 


(5) 


has been proposed as a direct measurement of the factorization breakdown effect and to ex¬ 
plore the px-dependent event plane angle fluctuations in the context of hydrodynamics ll23l . 
Here, the V n & coefficients are calculated by pairing particles within the same px interval (de¬ 
nominator) or from different px intervals (numerator). If the factorization relation of Eq. (|3j 
holds, this ratio is expected to be unity However, with the presence of a px-dependent event 
plane angle, it can be shown that the factorization ratio, r n , is equivalent to 


r n{Pj, Pt) 


(■ V n (pj)v n (p\) COS {n [Y„(Px) - fn(Pj)] }) 

\/( v UPt))(vI(Pt)) 


( 6 ) 


where Y n {Pj) arid Y n (pj) represent the event plane angles determined using particles from p'j 
and p'j intervals, respectively 11231 [24l. and () denotes averaging over all the events. As one can 
see from Eq. j6|, r„ is in general less than unity in the presence of the pj-dependent event plane 
angle Y„ fluctuations. 


4.2 Results for PbPb data 

The first measurement of pj-dependent factorization breakdown in PbPb collisions was pre¬ 
sented in Ref. fl3l . Our analysis is expanded to cover a much wider centrality range from 0% 
to 50%, and also includes a systematic comparison to hydrodynamic models. The values of 
r i(Pj/ Pt) and P3 (Pt/Pt) in PbPb collisions at v /s^ = 2.76 TeV are presented as a function of 
Pt — Pt ' n Figs. [I] and [ 2 } for several p‘j ranges in seven different centrality classes from 0-0.2% 
to 40-50%. The average px values within each p‘j and Pj range are used in order to calculate 
the difference between p‘j and p\. By construction, the r n value for the highest analyzed p\ 
range, where both particles are selected from the same pj interval, is equal to one. Only re¬ 
sults for p'j > Pj are presented, with a maximal p‘j value of 3 GeV/c, a kinematic regime where 
the hydrodynamic flow effect is believed to be dominant. The error bars correspond to statis¬ 
tical uncertainties, while systematic uncertainties are found to be negligible for the r n results 
(mainly because systematic uncertainties of V n & are typically on the order of a few percent, and 
ratios of V n ,\ are taken to form r n in this paper, where systematic uncertainties mostly cancel) 
and thus are not shown in any of the figures. 

A clear deviation from unity of the r 2 value (Fig. [TJ is observed for the highest pj ranges in 
very central PbPb collisions. For each centrality class, the effect becomes more pronounced 
with an increase of p‘j and also the difference between p l j and Pj values. This trend is expected 
as event-by-event initial-state geometry fluctuations play a more dominant role as the colli¬ 
sions become more central. The factorization breakdown effect reaches 20% in the ultra-central 
0-0.2% events for the greatest difference between pj and pj. For more peripheral centrality 
classes, the maximum effect is a few percent. Calculations using viscous hydrodynamics 6241 
are performed in all centrality classes, and shown as the curves in Fig. |T] To focus on the effect 
of initial-state fluctuations, the p/s value is fixed at 0.12. Two different models of initial condi¬ 
tions, MC-Glauber [46.471 and MC-KharzeevFevinNardi (MC-KFN; motivated by the concept 
of gluon saturation) |48|. are compared to data. The qualitative trend of the data is consistent 
with hydrodynamic calculations. However, quantitatively, neither of the two models can de¬ 
scribe all the data. The MC-Glauber model matches better the data for central collisions, while 
MC-KFN model appears to describe the data in the peripheral centrality range. 
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1.0 < p* < 1.5 GeV/c 1.5 < p* < 2.0 GeV/c 2.0 < p* < 2.5 GeV/c 2.5 < p* < 3.0 GeV/c 



Figure V. (Color online) The /^'dependent factorization ratio, ri, as a function of Pj — Pj in bins 
of pj for different centrality ranges of PbPb collisions at v /s NX = 2.76 TeV. The curves show 
the calculations from a viscous hydrodynamic model |24] using MC-Glauber and MC-KLN 
initial condition models, and an p/s value of 0.12. Each row represents a different centrality 
range, while each column corresponds to a different p‘j range. The horizontal solid lines denote 
the f 2 value of unity The error bars correspond to statistical uncertainties, while systematic 
uncertainties are negligible for the r„ results, and thus are not shown. 


For the third-order harmonics (n = 3), the effect of factorization breakdown is significantly 
smaller than for the second-order harmonics. Only a weak centrality dependence of r$ is seen 
in Fig. [ 2 ] The biggest deviation of r 3 from unity is about 5% at large values of Pj — p^. (i.e., > 
1 GeV/c). Again, the qualitative features of the data are described by the hydrodynamic model, 
although the effects are over-estimated for peripheral collisions by the model. Calculations of 
7*3 using two different initial-state models yield similar results, with MC-KFN model showing 
a slightly stronger centrality dependence. 

To understand better how the effects of factorization breakdown and pi-dependent event plane 
fluctuations are influenced by the initial-state conditions and the value of rj / s in hydrody- 
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4 Transverse momentum dependence of factorization breakdown 
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Figure 2: (Color online) Similar distributions as shown in Fig. [lj but for the factorization ratio 


namic models, a detailed comparison of measured r 2 values in 0-0.2% centrality PbPb colli¬ 
sions (where the effect is most evident) to hydrodynamic calculations is shown in Fig. [3] For 
this comparison, calculations with MC-Glauber and MC-KLN initial conditions are each per¬ 
formed for three different ?/ / s values and compared to data. For each initial-state model, the )'2 
values are found to be largely insensitive to different values of ?//s. This is because, in defining 
r„(pj, Pj), the magnitudes of anisotropy harmonics, which have a much greater sensitivity to 
i) / s, are mostly canceled. Fluctuations of the event plane angle in pj are mainly driven by the 
non-smooth local fluctuations in the initial energy density distribution. This comparison shows 
that the use of r n data can provide new constraints on the detailed modeling of the initial-state 
condition and the fluctuations of the medium created in heavy ion collisions, which is indepen¬ 
dent of the i] /s value. The better constraints on the initial-state conditions found using the r n 
data will, in turn, improve the uncertainties of determining the medium's transport properties 
(e.g., ;//s) using other experimental observables (e.g., the v n magnitude, which is sensitive to 
both the initial state and f]/s). 
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1.0 < p° < 1.5 GeV/c 1.5 < p“ < 2.0 GeV/c 2.0 < p“ < 2.5 GeV/c 2.5 < p‘ < 3.0 GeV/c 



Figure 3: (Color online) Factorization ratio, r^, as a function of p‘j — Pj in bins of Pj for 0-0.2% 
centrality PbPb data at = 2.76 TeV compared to viscous hydrodynamic calculations l24l 
using MC-Glauber and MC-KLN initial condition models, and three different values of / 7 /s. 
The horizontal solid lines denote the r 2 value of unity. The error bars correspond to statistical 
uncertainties, while systematic uncertainties are negligible for the r n results and thus are not 
shown. 
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Figure 4: (Color online) The pi-dependent factorization ratio, ri, as a function of Pj — Pj in bins 
of p‘j for four N°® me ranges in 5.02 TeV pPb collisions. The curves show the predictions from 
hydrodynamic calculations for pPb collisions of Ref. 11251 . The horizontal solid lines denote 
the r 2 value of unity. The error bars correspond to statistical uncertainties, while systematic 
uncertainties are negligible for the r n results and thus are not shown. 


To gain insights into the origin of long-range correlations observed in high-multiplicity pPb 
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4 Transverse momentum dependence of factorization breakdown 


1.0 < < 1.5 GeV/c 1.5< p‘< 2.0 GeV/c 2.0 < p% 2.5 GeV/c 2.5 < p" < 3.0 GeV/c 



Figure 5: (Color online) Similar distributions as shown in Fig. [Ij but for the factorization ratio 

n 3 - 

collisions, the measurement of and rj is also performed for pPb data at ^/s NN = 5.02 TeV 
for four different high-multiplicity ranges. The results are shown in Figs. [4] and [5j in the same 
format as those for PbPb collisions, for four pj ranges (of increasing pj from left to right panels) 
as a function of pj — p\. 

Breakdown of factorization is observed in the ip results of pPb collisions for all multiplicity 
ranges investigated in this paper. Similar to PbPb collisions, for any multiplicity range, the 
effect gets larger with an increase in the difference between pj and Pj values. However, the 
observed factorization breakdown reaches only up to 2-3% for the largest value of p‘j — p\ at 
2.5 < pj < 3.0 GeV/c. This is significantly smaller than that seen in central PbPb collisions. 
Little multiplicity dependence of r 2 is observed in pPb collisions. Comparison of the CMS data 
to hydrodynamic predictions for pPb collisions in Ref. Il25l is also shown. In this hydrody¬ 
namic calculation, a modified MC-Glauber initial-state model is employed for pPb collisions 
where the contributing entropy density of each participating nucleon in the transverse plane 
is distributed according to a 2D Gaussian distribution. The width of the transverse Gaussian 
function can be chosen to vary the transverse granularity of fluctuations, to which the r n values 
are found to be most sensitive. The r 2 data are better described by calculations with a width 
parameter of 0.4 fm (curves in Fig. |4|, while a width of 0.8 fm gives an r n value of nearly unity 
(not shown) and thus underestimates the effect observed in the data. For both cases, the cal¬ 
culations are found to be insensitive to different p /s values, consistent with the hydrodynamic 
calculations used for more central PbPb collisions presented earlier. 

Results of r 3 are shown in Fig.[5j presented in the same format as for r 2 - Within current statistical 









































4.4 Comparison of pPb and PbPb data 
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precision, no evident breakdown of factorization is found in very-high-multiplicity pPb events 
(185 < N°^ me < 260), while the r 3 value exceeds unity for much lower-multiplicity pPb events 
at high pj, particularly for 120 < N°® me < 150. This is a clear indication of significant nonflow 
effects as the event multiplicity decreases, because the r n values predicted by hydrodynamic 
models with px-dependent event plane angle fluctuations would always be equal to or less 
than one, according to Eq. (| 6 |. One obvious possibility is back-to-back jet correlations, which 
would give a large negative contribution to V 3 A at high p‘j and pj values in low multiplicity 
events IflOl . This would lead to a significant reduction of the denominator of Eq. (j6]> and drives 
the 73 value up above unity Very little effect of factorization breakdown for n = 3 is predicted 
in Ref. l25l . which is consistent with the data except for the low-multiplicity ranges. 


4.4 Comparison of pPb and PbPb data 


l.0<p"<l.5 GeV/c 
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Figure 6: (Color online) The px-dependent factorization ratio, ri, as a function of p‘j — Pj in 
bins of p'j for two N°® me ranges of 5.02 TeV pPb and 2.76 TeV PbPb collisions. The horizontal 
solid lines denote the ri value of unity. The error bars correspond to statistical uncertainties, 
while systematic uncertainties are negligible for the r n results and thus are not shown. 

Figure [6] shows a comparison between 5.02 TeV pPb and 2.76 TeV peripheral PbPb collisions 
over the same multiplicity ranges. Because of the statistical limitation of the PbPb data, the 
multiplicity ranges used in Figs. [I] and [ 5 ] for pPb data are combined into two N°® me classes, 
100 < < 185 (top) and 185 < N°® me < 260 (bottom). At a similar N°® me range, 

the magnitudes of factorization breakdown in pPb and PbPb collisions depart from unity by 
less than 8%, with slightly smaller deviations for pPb data, although the statistical precision 
is limited. For both high-multiplicity pPb and peripheral PbPb collisions, the observed effect 
is significantly smaller than that for 0-0.2% centrality ultra-central PbPb collisions (up to 20% 
away from unity). The similar behavior (e.g., px-dependence) of factorization data in pPb to 
that in PbPb collisions may provide new insight into the possible hydrodynamic flow origin of 
long-range two-particle correlations in the pPb system, particularly in providing new informa¬ 
tion on the nature of initial-state fluctuations in a much smaller volume. 


To study directly the multiplicity dependence of the effect in PbPb and pPb collisions, the r 2 
and 7*3 results for 2.5 < Px < 3.0 GeV/c and 0.3 < p\ < 0.5 GeV/c (where the difference between 
p l j and p ! j is the greatest, pj — Pt ~ 2GeV/c) are shown in Fig. as a function of event mul¬ 
tiplicity in pPb and PbPb collisions. Here, the number of tracks is still counted with |?/| < 2.4 
and px > 0.4 GeV/c but corrected for the detector inefficiency, since a different track recon- 



















12 


4 Transverse momentum dependence of factorization breakdown 
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Figure 7: (Color online) The px-dependent factorization ratios, 72 and r 3 , as a function of event 
multiplicity in pPb and PbPb collisions. The curves show the calculations for PbPb collisions 
from viscous hydrodynamics in Ref. Il24l with MC-Glauber and MC-KLN initial condition mod¬ 
els and 77 /s = 0.12, and also hydrodynamic predictions for PbPb and pPb data in Ref. H25l . The 
horizontal solid lines denote the i '2 (top) and 73 (bottom) value of unity. The error bars corre¬ 
spond to statistical uncertainties, while systematic uncertainties are negligible for the r n results 
and thus are not shown. 

struction algorithm is used for the pPb and central PbPb data. Additionally, at the top of the 
figure, a centrality axis is shown which is applicable only to PbPb collisions. The breakdown of 
factorization for r 2 in PbPb events increases dramatically as the collisions become more central 
than 0-5%, while the effect in 73 remains at the 2-3% level, largely independent of centrality. 
For more peripheral PbPb events from 20% to 80% centrality, the deviation of 72 from unity 
increases slightly from about 2 to 5%. Calculations using a hydrodynamic model in PbPb colli¬ 
sions ||24| with MC-Glauber and MC-KLN initial conditions and rj/s = 0.12 are also shown as 
dotted and dash-dotted curves, respectively, as a function of centrality. As pointed out earlier, 
neither of the two calculations can describe the data quantitatively over the entire centrality 
range, although the qualitative trend is reproduced. The 72 values for pPb show little multi¬ 
plicity dependence, consistent with hydrodynamic predictions in Ref. Il25l . The 73 values for 
pPb go significantly above unity at lower multiplicities, because of the onset of nonflow corre¬ 
lations. The discrepancy in the hydrodynamic calculations of 72 for peripheral PbPb collisions 
between Ref. Il24)l and Ref. l25l may be related to differences in some model parameters (e.g., 
transverse size of the nucleon). This should be investigated in the future. 

Although the factorization results presented in this paper suggest a breakdown of the assump- 
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tion commonly applied in studying collective flow using two-particle correlations (Eq. ([3])), 
previous v n measurements from the two-particle method still remain valid. However, they 
should be more precisely interpreted as the v n values obtained with respect to an averaged 
event plane using particles from a given kinematic regime (usually over a wide pj range). The 
studies in this paper also point out the importance of applying the same conditions for theoret¬ 
ical calculations when comparing to the experimental data. 


5 Pseudorapidity dependence of factorization breakdown 

5.1 Analysis technique 

In principle, the //-dependent factorization breakdown and event plane angle fluctuations can 
be examined using a similar formalism to Eq. ^ by replacing p‘j and p\ with if and if. How¬ 
ever, the main issue with this approach is that the requirement of |A//| > 2 for removing 
short-range two-particle correlations cannot be fulfilled anymore as the denominator of the 
factorization ratio takes the V„A{p a , if) components, where //" ~ if. The correlation signal 
from collective flow is strongly contaminated by short-range jet-like correlations. To avoid this 
problem, an alternative observable is developed for the study of //-dependent factorization, by 
taking advantage of the wide // coverage of the CMS tracker and HF calorimeters. 

The //-dependent factorization ratio, r n (tj a , if), is defined as 


r n {v a > f l b ) 


VnA{-lf,lf) 
VnAff^f) ' 


(7) 


where V n ^(i] a , if) is calculated in the same way as Eq. Q but for pairs of particles taken from 
varied if and if regions in fixed p‘j and p\ ranges. Here, particle a is chosen from charged 
tracks with 0.3 < pj < 3.0GeV/c and ;/"| < 2.4, while particle b is selected from the HF cal¬ 
orimeter towers with the energy exceeding 1 GeV (with a total coverage of 2.9 < |//| < 5.2) 
without any explicit transverse energy ( Ej) threshold for each tower. With this approach, the // 
values of both particles from a pair can be varied over a wide range, while it is possible to en¬ 
sure a large // gap by combining detector components covering central and forward // regions. 
As illustrated by a schematic in Fig. [8j for 4.4 < if < 5.0 from the HF calorimeters, a minimum 
// gap of 2 units between a calorimeter tower and any charged particle from the silicon tracker 
is guaranteed. Away-side back-to-back jet correlations could still be present but they are shown 
to have a negligible contribution at low pj because of very high multiplicities l22l . especially 
in central PbPb collisions. To account for any occupancy effect of the HF detectors due to large 
granularities in // and cp, each tower is weighted by its Ej value when calculating the average 
in Eq. j4j. For consistency, each track is also weighted by its pj value. The finite azimuthal 
resolution of the HF towers (0.349 radians) has negligible effects to the V n ,\ calculation, which 
takes an Ej-weighted average of 36 tower segments over a 2rc coverage. 

If, for each event, the event plane angle, Y„, does vary for particles produced at different // 
regions, the following relation for r n (if, if) can be derived. 


, a jjb) = (Vn(->f)Vn(lf)cOs{lI [Y n (-lf)-Y n (lJ b )]}) 

(v„(if)v n (iJ b )cos{n [Y„(/f) -Y n (tf)\j) 


In symmetric collision systems like PbPb, v n harmonics from symmetric positive (Vn(if )) and 
negative (v n (—?/")) // regions are identical after averaging over all events. Therefore, Eq. (js]) can 
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5 Pseudorapidity dependence of factorization breakdown 


rnOf.Tf) = 


V nA (Tf,T1 b ) 





V nA (r, a ,r, b ) 


Tracker 



-5.2 


-3.0 -2.4 


2.4 3.0 


11 


5.2 


Figure 8: (Color online) A schematic illustrating the acceptance coverage of the CMS tracker 
and HF calorimeters, and the procedure for deriving the //-dependent factorization ratio, 

r n (ri a ,i] b ). 


be approximated by 


r n {rj a ,if ) 


(cos [n (Y„(-// fl ) - Y„(// b ))]; 


(9) 


(cos [n (Y,,(//") — Y„(// fc ))]) 

Here, the approximation is due to the fact that the flow magnitude (v n ) and the orientation 
angle (Y ;1 ) are inside the same averaging over all the events in the numerator of Eq. (j8|. As a 
result, r n {rj a , if) represents a measurement of relative event plane angle fluctuations in // for 
planes separated by |// fl + t] b \ and \rj a — i) b \. Similar to r n (pj, p\), r n (p a , if) is equal to unity if 
the factorization holds but factorization breaks down in general in the presence of event plane 
fluctuations in //. 

For an asymmetric collision system like pPb, v n (t] a ) and "o n ( — rf ) are not identical in general, 
and thus //-dependent event plane fluctuation effect cannot be isolated in Eq. However, by 
taking a product of r n (rj a , if) and r n ( — ij a , —rj b ), the v„ terms can be removed, 

Jr 77a „h )r 7 ~a „^^,/ ( cos [”( Y »(-V fl )- Y »(V fc ))]) (cos [n (Y„(>/' ; ) — Y„(—?/ b ))]) 
V""'' ” V ( cos [n (Y fl (i/“) - Y„(// fc ))]) (cos [n (Y„(-// fl ) - Y„(-// b ))])' 

( 10 ) 

In this way, the //-dependent event plane angle fluctuations in pPb collisions can also be stud¬ 
ied. 

5.2 Results for PbPb data 

The results of //-dependent factorization ratios, f 2 , r$ and r 4 in PbPb collisions at = 

2.76 TeV are shown in Figs. [9]- 11 as a function of //" for eight different centrality classes from 
0-0.2% to 50-60% (except for /q, for which only three centrality classes are shown because of 
statistical limitations). The zq {p n ,p b ) values are calculated in //" bins of 0.3 units, and the if 
value at the center of each bin is used in the plots. Data obtained with calorimeter tower // 
ranges 3.0 < i] b < 4.0 and 4.4 < rj h < 5.0 are both presented. Since PbPb is a symmetric sys¬ 
tem, the V n (rj a , if) and V n (—rj a , —if) coefficients are combined before calculating the r n ratios 
in order to achieve the optimal statistical precision. Charged tracks within 0.3 < pj < 3.0 GeV/c 
and all calorimeter towers (E > 1 GeV) are used. When //" = 0, the r n value is equal to unity 
by construction since both the numerator and denominator of r n have the same // gap between 
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Figure 9: (Color online) The //-dependent factorization ratio, 7 * 2 , as a function of //" for 3.0 < 
< 4.0 and 4.4 < if < 5.0, averaged over 0.3 < pj < 3.0GeV/c, in eight centrality classes of 
PbPb collisions at ^/snn = 2.76 TeV. The curves correspond to fits to the data for 4.4 < if < 
5.0 given by Eq. 0 - The horizontal solid lines denote the r 2 value of unity The error bars 
correspond to statistical uncertainties, while systematic uncertainties are negligible for the r„ 
results, and thus are not shown. 
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particle a and b, as indicated in Eq. (j9ji. As rj a increases, a significant decrease of r n below unity 
is observed, which may suggest the presence of //-dependent event plane angle fluctuations. 

The r 2 values for 4.4 < i] 1 ’ < 5.0 are found to decrease with //" approximately linearly for most 
of the centrality classes up to a few percent deviation below unity at rj a ~ 2.4. This behavior 
is slightly different for the most central 0 - 0 . 2 % events, where the decrease of V 2 becomes more 
significant at if 1 ~ 1. For 3.0 < //’ < 4.0, the r 2 value exhibits a much stronger factorization 
breakdown effect for an if > 1. This can be understood as the effect of short-range jet-like 
correlations when the // gap between two particles is less than 2 , which increases the denomi¬ 
nator of Eq. |7]). However, for if < 1, the V 2 results are found to be consistent with each other, 
independent of rj b (except for 0-0.2% centrality). This demonstrates that contributions of short- 
range jet-like correlations are almost completely suppressed if the requirement of \Arj\ > 2 to 
both numerator and denominator of r n (rj a , if) is imposed. 


The effect of //-dependent factorization breakdown is much stronger for higher-order harmon¬ 
ics, ro, and /q, shown in Figs. 10 and 11 For /q, this trend is opposite to what is observed for 
the pi-dependent factorization ratio. For all centrality ranges (including 0-0.2%), an approx¬ 
imate linear dependence of r 3 and r 4 is seen. Results from the two different if ranges agree 
over most of the if range within statistical uncertainties. This might suggest that short-range 
jet-like correlations have much smaller effects on higher-order harmonics. 

As observed in Figs. [9 - 11 the r n ( yj a , if) values are independent of if, for if ranges where con¬ 


tributions of only long-range (|A//| > 2) correlations are included. To quantify the dependence 
of r„ values on if, a simple empirical parameterization is introduced: 


n (f¥ n (rj a )-Y n (r, b )) 


cos 


= e -F V n\v“-V b \ / 


( 11 ) 
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5 Pseudorapidity dependence of factorization breakdown 



Figure 10: (Color online) Similar distributions as shown in Fig. [9j but for the factorization ratio 

h3- 



Figure 11: (Color online) Similar distributions as shown in Fig. [ 9 J but for the factorization ratio 
jq in fewer centrality ranges. 


which is based on the assumption that relative fluctuations between two event plane angles 
depend only on their pseudorapidity difference. At small At] values, the exponential function 
form can be approximated by a linear function of At], consistent with the observation in the 
data. By plugging Eq. ( fTTj l into Eq. the r„ can be expressed as 


r„{if,r] h 


-2F„y 


( 12 ) 


11 


According to Eqs. (Ill 


which is independent of t] b , consistent with the results in Figs. [9- 

and (12 > the r n (//", ?/' ) also corresponds to a measurement of event plane fluctuations between 

Y„()f) and Y n (-t] a ), 


r n (v a >V b ) ~ ( cos [n (Y„(-7 fl ) - Y„(?/ fl ))]). 


(13) 


The r-i data for 4.4 < tf < 5.0 are well fit with a functional form given by Eq. (12 1 for most 
centrality classes (x 2 /(degree of freedom) ~1), except for 0-0.2% centrality, where the rj value 































5.3 Results for pPb data 
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deviates from unity much faster as if increases. Note that the parameter, F//, is purely empir¬ 
ical, without any clear physical meaning at present. It is introduced mainly for quantitatively 
evaluating the centrality evolution of factorization breakdown effect, as will be discussed later 
in Section l5~4l 

5.3 Results for pPb data 



Figure 12: (Color online) The square root of the product of factorization ratios, 
\/r 2 (t] a ,i] b )r 2 (—t] a ,—i] b ), as a function of tj a for 3.0 < rf < 4.0 and 4.4 < if < 5.0, averaged 
over 0.3 < pj < 3.0 GeV/c, in four multiplicity classes of pPb collisions at ^ = 5.02 TeV. 
The curves correspond to the fits to the data for 4.4 < rf < 5.0 using Eq. (n2k. The horizontal 
solid lines denote the r 2 value of unity The error bars correspond to statistical uncertainties, 
while systematic uncertainties are negligible for the r 2 results, and thus are not shown. 


Studies of //-dependent factorization breakdown of two-particle correlations are also performed 



thus does not only reflect the effect of event plane angle fluctuations. Therefore, the results are 
presented as the square root of the product of r n (i] a , if) and r„(— rj a , —if), which is designed 
to remove the sensitivity to the magnitude of v n (see Eq. (To} for details). Similar to those in 
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6 Summary 


PbPb collisions, two different ;/ ranges of HF towers, 3.0 < t] b < 4.0 and 4.4 < t] h < 5.0, are 
compared. 

A significant breakdown of factorization in ij is also observed in pPb collisions as rj a increases. 
Similar to the PbPb results, the factorization breakdown is approximately independent of ij b for 
i) a < 1 for all multiplicity ranges but shows a much larger deviation from unity for 3.0 < tj b < 
4.0 as if increases beyond one unit because of short-range correlations. The fits to the data for 
4.4 < if < 5.0 using Eq. (12 1 are also shown; the data are well-described over the accessible i] n 


range. It should be noted that the assumption made in Eq. (Ill is purely an empirical param¬ 
eterization for quantifying the behavior of the data. Since pPb collisions are asymmetric, this 
assumption could be invalid. More detailed investigations on how r n depends on if 1 and if in 
the proton- and lead-going directions, respectively, are needed in future work. 

5.4 Comparison of pPb and PbPb data 


The extracted F„ parameters are plotted as a function of event multiplicity in Fig. 13 in pPb 


collisions for n = 2 and PbPb collisions for n = 2-4. The Fj value reaches its minimum around 
midcentral (~20%) PbPb events, and increases significantly for more peripheral PbPb events 
and also for pPb events, where the relative fluctuations of u>i are larger fT2f. Toward the most 
central PbPb events, the Fj value also shows a tendency to increase slightly, although the r n 
data for 0-0.2% centrality are not well described by Eq. ( [12) . At a similar multiplicity, magni¬ 
tudes of the Fj parameter in pPb are significantly larger than those in PbPb, and decrease with 
increasing event multiplicity. In PbPb collisions, a much stronger if -dependent factorization 
breakdown is seen for higher-order harmonics than for the second order, as shown by the F| 
and F^ parameters. There is little centrality dependence for n = 3, except for the most cen¬ 
tral 0-20% PbPb collisions. Within current statistical uncertainties, no centrality dependence is 
observed for n = 4. 


6 Summary 

Factorization of azimuthal two-particle correlations into single-particle anisotropies has been 
studied as a function of transverse momentum and pseudorapidity of each particle from a pair, 
in PbPb collisions at ^/s NN = 2.76 TeV and pPb collisions at ^ = 5.02 TeV, over a wide multi¬ 
plicity range. The factorization assumption is found to be broken as a function of both pj and ;/. 
The effect of pj dependent factorization breakdown for the second-order Fourier harmonic is 
found to increase with the difference in pj between the two particles. The factorization break¬ 
down reaches 20% for the most central PbPb collisions, while it decreases rapidly for more 
peripheral collisions. The effect is significantly smaller (2-3%) in high-multiplicity pPb colli¬ 
sions. In both PbPb and pPb samples over the full centrality or multiplicity range, little effect 
is observed for the third order harmonic. For the // dependence, the observed factorization 
breakdown shows an approximately linear increase with the ?/ gap between two particles for 
all centrality and multiplicity classes in PbPb and pPb collisions. The effect is weakest for mid¬ 
central PbPb events, but becomes larger for more central or peripheral PbPb collisions, and also 
for very high-multiplicity pPb collisions. Moreover, a much stronger //-dependent effect is seen 
for the third- and fourth-order harmonics than the second-order harmonics in PbPb collisions. 
This relation between the second and third order is opposite to that seen in the pj-dependent 
factorization studies. The observed factorization breakdown presented here does not invalidate 
previous v„ measurements. Instead, the previous values should be reinterpreted as measuring 
anisotropies with respect to the event plane averaged over a given kinematic region. Further¬ 
more, it is important to compare data and theoretical calculations following exactly the same 
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PbPb centrality(%) 



Figure 13: (Color online) The Fjj parameter as defined in Eq. (12 ' as a function of event multi¬ 
plicity in PbPb collisions at v /s^" = 2.76 TeV for n = 2-4 and pPb collisions at ^ = 5.02 TeV 
for n = 2. The error bars correspond to statistical uncertainties, while systematic uncertainties 
are negligible for the r n results, and thus are not shown. 


procedure. 

The factorization data have been compared to hydrodynamic calculations with fluctuating 
initial-state conditions. The pj-dependent factorization data are qualitatively described by vis¬ 
cous hydrodynamic models, which are shown to be largely insensitive to the value of shear 
viscosity to entropy density ratio of the medium. This observation offers great promise for 
using the factorization data to disentangle contributions of the initial-state conditions and the 
medium's transport properties to the observed collective flow phenomena in the final state. 
The new studies of / 7 -dependent factorization breakdown give an indication of initial-state 
fluctuations along the longitudinal direction. This will provide new insights into the longi¬ 
tudinal dynamics of relativistic heavy ion collisions, and help improve the three-dimensional 
modeling of the evolution of the strongly-coupled quark gluon medium. 
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